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Abstract 


Prevalence of Middle East respiratory syndrome coronavirus (MERS-CoV) was determined in 45 dromedary camels 
from two geographically separated herds in Jordan. Virus shedding was only detected in swabs obtained from the 
respiratory tract and primarily observed in camels younger than 3 years. MERS-CoV seroprevalence increased with 
age of camels. Bovine and sheep sera were seronegative. Phylogenetic analysis of partial S2 clustered the Jordanian 


MERS-CoV strains with contemporary MERS-CoV strains associated with nosocomial outbreaks. 
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Introduction 


S INCE THE IDENTIFICATION of Middle East respiratory 
syndrome coronavirus (MERS-CoV) as the causative 
agent of a fatal case of respiratory tract disease in the Kingdom 
of Saudi Arabia (KSA) in 2012, the virus has caused >1700 
laboratory-confirmed cases of disease, including >600 fatal 
cases (WHO 2016). The first known MERS-CoV outbreak in 
humans, diagnosed retrospectively, occurred in Jordan in 2012 
(Hijawi et al. 2013). The dromedary camel has been identified 
as the primary reservoir of MERS-CoV, and direct evidence for 
zoonotic transmission from camels has been reported in KSA 
and Qatar (Haagmans et al. 2014, Memish et al. 2014). A 
limited number of seropositive camels in Jordan have been 
described previously (Reusken et al. 2013), but no MERS-CoV 
has been detected in camels in Jordan. 


Materials and Methods 


During May 2016, we collected swabs (from nasal, uro- 
genital, and rectal areas) and blood samples from camels at 


two locations in Jordan. Two camel herds were identified to 
study (1) a traditional bedouin camel herd, in which camels 
are allowed to graze and browse freely (Azraq, Zarqa) and (2) 
a more conventional mixed farm setting, where camels were 
kept in pens on one farm (Ramtha) (Fig. 1A). In Ramtha, 
blood samples from adult sheep and cattle were additionally 
collected. Swab samples were collected in virus transport 
medium. RNA was extracted from samples using the 
QiaAmp Viral RNA kit (Qiagen). Five microliters of RNA 
was used in a one-step real-time RT-PCR UpE assay for 
MERS-CoV using the Rotor-Gene™ probe kit (Qiagen). 
Positive samples (cycle threshold [Ct] < 37) were tested using 
the ORFIA assay (Corman et al. 2012), and samples were 
excluded from further analysis when ORFIA testing was 
negative. cDNA was synthesized using random hexamer and 
used to PCR amplify the MERS-CoV spike S2 domain (nu- 
cleotides 23781-24395 of HCoV-EMC/2012) as described 
previously (Smits et al. 2015). Sequences were assembled 
on SeqMan Pro and analyzed on MegAlign (DNASTAR). 
Phylogenetic trees of the S2 domain were generated using 
Mega 6.0.6 with the maximum likelihood statistical method 
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MERS-CoV prevalence in dromedary camels in Jordan. (A) Locations of dromedary camels (north: Ramtha; east: 


Azraq). (B) Percentage of MERS-CoV RNA shedding dromedary camels as detected by UpE and ORFIA qRT-PCR assay 
in nasal swab, stratified by age. (C) Percentage of MERS-CoV S1-specific seropositive dromedary camels, stratified by age. 
(D) ELISA ratio of seropositivity of dromedary camels, stratified by age. The ELISA ratio was calculated by dividing the 
OD of each serum sample by a constant positive sample on the ELISA plate. ELISA, enzyme-linked immunosorbent assay; 
Middle East respiratory syndrome coronavirus (MERS-CoV), OD, optical density. 


based on the GTR+G+I model with 1000 bootstraps repli- 
cates. Sera were analyzed by MERS-CoV spike protein (S) 
enzyme-linked immunosorbent assay (ELISA); Maxisorp 
(Nunc) plates were coated overnight with S1 protein (Sino 
Biological) and blocked with 1% milk. Sera (400x dilution) 
were added to the plate in duplicate. MERS-CoV S1- 
specific antibodies were detected using anti-llama (Agri- 
sera), anti-bovine, or anti-sheep (Jackson) IgG (H&L) 
HRP-conjugated antibody on development with peroxidase- 
substrate reagent (KPL); optical density (OD) was measured 
at 405 nm. The threshold of positivity was mean OD + 10x 
standard deviation of negative sera obtained from camels 
raised in captivity in the United States. 


Results 


Twenty-three camels sampled from the Bedouin herd in 
Azraq ranged in age from 4 months to 8 years, whereas 22 
camels sampled at the farm in Ramtha ranged in age from 4 
months to 3 years (Table 1). Ten sheep and five cows were 
sampled. Nasal discharge was observed at the time of sam- 
pling in some of the MERS-CoV-positive camels <1 year old 
but not in older camels. 

RNA from 42/45 camels was tested positive for the pres- 
ence of MERS-CoV nucleic acid. MERS-CoV RNA was 


solely detected in nasal swabs. Importantly, urogenital and 
rectal swab samples were negative. One nasal swab origi- 
nating from camel 40 was excluded from analysis; while the 
UpE assay resulted in a Ct value of <37, the ORFIA assay 
was negative. An unpaired two-tailed Student’s f-test com- 
paring the age of viral RNA-positive versus RNA-negative 
camels was significant (p=0.0311). Only 1/7 camels older 
than 3 years was positive for viral RNA, whereas 11/18 
camels <1 year old, 4/4 camels 1-2 years old, and 12/13 
camels 2-3 years old were positive for viral RNA in nasal 
swabs (Fig. 1B). MERS-CoV-specific antibodies were ob- 
served in the majority of animals: 78% of camels <1 year old, 
69% of camels 1—2 years old, and 100% of camels >2 years 
old were seropositive (Fig. 1C). No MERS-CoV S1-specific 
antibodies were found in sheep or bovine serum samples. 
To better assess the potential significance of ELISA val- 
ues as a correlate of MERS-CoV susceptibility, we calcu- 
lated the ratio of the ELISA value of each sample to that 
obtained from an included camel reference serum on each 
plate in the assay. By this approach, we show that camels <3 
years old exhibited an average ELISA ratio of 0.5-0.65, 
whereas older animals exhibited an average ELISA ratio of 
1.63. This difference was significant as tested via a two- 
tailed Mann—Whitney test (Fig. 1D). Spike S2 partial domain 
sequences were obtained from 16/28 samples (accession 
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TABLE |. AGE, LOCATION, AND ASSAY RESULTS PER DROMEDARY CAMEL 


TCID50eq/swab 
Camel Age Herd Virus positive UpE ORFIA Seropositive ELISA ratio 
40 2M Ramtha 0 1 1.15 
41 2M Ramtha 0 1 0.50 
42 2M Ramtha 1 135 156 0 0.11 
07 4M Azraq 1 12 37 0 0.18 
13 4M Azraq 1 15 93 0 0.19 
33 4M Ramtha 1 252 196 1 0.17 
36 4M Ramtha 0 1 0.79 
09 5M Azraq 1 200 728 1 0.43 
12 5M Azraq 1 21 35 0 0.21 
16 5M Azraq 0 1 0.35 
35 5M Ramtha 1 92 109 1 0.29 
06 6M Azraq 1 81 264 1 0.30 
17 6M Azraq 1 188 174 1 1.27 
24 6M Ramtha 0 1 1.43 
34 6M Ramtha 0 1 0.43 
37 6M Ramtha 1 2729 1256 1 0.30 
10 7™M Azraq 0 1 0.58 
11 7™M Azraq 1 70 209 1 1.07 
25 12M Ramtha 1 9575 4420 0 0.13 
26 13M Ramtha 1 1569 2643 0 0.09 
30 14M Ramtha 1 5139 6054 0 0.13 
32 14M Ramtha 1 5907 2760 0 0.09 
18 15M Azraq 1 56 802 1 1.46 
08 18M Azraq 1 76 86 1 1.25 
14 18M Azraq 1 128 371 1 1.05 
15 18M Azraq 1 23 119 1 1.35 
27 18M Ramtha 1 1067 1074 1 0.16 
28 18M Ramtha 1 507 876 1 0.15 
29 18M Ramtha 1 168 368 1 0.63 
31 18M Ramtha 1 491 230 1 0.21 
21 19M Azraq 0 1 1.79 
38 30M Ramtha 1 280 452 1 1.10 
39 30M Ramtha 1 152 155 1 0.54 
44 2Y Ramtha 1 112 157 1 0.57 
45 2Y Ramtha 1 76 219 1 0.15 
01 3Y Azraq 0 1 1.76 
20 3Y Azraq 0 1 1.74 
43 3Y Ramtha 1 87 179 1 0.82 
02 4Y Azraq 0 1 1.83 
04 SY Azraq 0 1 1.85 
22 SY Azraq 0 1 1.89 
03 6Y Azraq ND 1 1.76 
19 6Y Azraq 0 1 1.23 
05 8Y Azraq ND 1 1.81 
23 8Y Azraq ND 1 1.67 


ELISA ratio was defined by dividing the OD of each serum sample by a constant positive sample on the ELISA plate. All viral RNA was 


detected in respiratory tract swab samples, none in urogenital or fecal swab samples. Positivity is marked as ‘‘1,” 


negativity as “0.” 


ELISA, enzyme-linked immunosorbent assay; M, month; OD, optical density; Y, year. 


numbers KX443663-KX443678). We were unable to obtain 
sequences of 12 positive samples with a lower viral load (Ct 
>33). We performed phylogenetic analysis with a selection 
of MERS-CoV S2 sequences using representatives of known 
clades, as performed previously (Smits et al. 2015). The 
phylogenetic analysis of the partial S2 sequences placed the 
novel Jordanian viruses within the B1 cluster representing 
contemporary camel and human MERS-CoVs. Three se- 
quences were identical to spike S2 sequences of human 
isolates of MERS-CoV obtained in Jordan in 2015 (Lamers 


et al. 2016). Twelve samples differed by two synonymous 
mutations (C23837T, T24074G), whereas the remaining 
samples contained a mixture of 24074T and 24074G com- 
bined with 23837T. Importantly, all novel sequences were 
found in cluster B1, containing the most recent MERS-CoV 
sequences originating from camel and human viruses 
(Fig. 2). Thus, the circulating MERS-CoV strains in the 
Jordan dromedary camel population are closely related to 
virus strains known to be capable of zoonotic transmission 
and to cause disease in the human population. 
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FIG. 2. Phylogenetic analysis of a partial spike $2 domain. 
A maximum likelihood tree based on the GTR+G+I model 
using 1000 bootstraps was generated from a spike S2 domain 
genome fragment corresponding to nucleotides 23781-24395 
of HCoV-EMC/2012. The newly identified MERS-CoV se- 
quences are depicted in bold, recent MERS-CoV sequences 
associated with an outbreak in Jordan in 2015 are depicted in 
bold. Bootstrap values of <50 were not shown. 
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Discussion 


This study confirms the circulation of MERS-CoV within the 
dromedary camel population of Jordan in line with MERS-CoV 
detection in camels throughout the Middle East (Reusken et al. 
2016). Importantly, MERS-CoV RNA could be detected in nasal 
swabs of seropositive dromedary camels. Antibodies against 
MERS-CoV were found at a young age in dromedary camels. 
The presence of maternal antibodies via the intake of colos- 
trum during the first 24h postparturition could play a role in 
the detection of MERS-CoV-specific antibodies in animals 
<6 months old and might not reflect actively acquired anti- 
bodies (Kamber et al. 2001, Meyer et al. 2016). 

We observed shedding of MERS-CoV in the presence of 
antibodies, which suggests either reinfection of seropositive 
animals or shedding of virus/viral RNA during early stages of 
seroconversion. Previously, a lack of correlation was ob- 
served between virus/viral RNA shedding and the presence of 
neutralizing antibodies; these data highlight the potential for 
reinfection of seropositive animals (Farag et al. 2015). The 
ELISA ratio might be a better predictor of MERS-CoV sus- 
ceptibility than seropositivity; 24/29 animals with an ELISA 
ratio of <1.1 were MERS-CoV positive, whereas only 4/16 
animals with an ELISA ratio of 21.1 were positive for viral 
RNA in nasal swabs. In contrast, 20/34 seropositive drome- 
dary camels were MERS-CoV viral RNA positive versus all 
(8/8) seronegative dromedary camels. This might indicate 
that sterile immunity is only reached at high levels of anti- 
body titers, in line with a previous study conducted in the 
United Arab Emirates (Meyer et al. 2016). 

The phylogenetic analysis of the partial S2 sequences 
placed the circulating viruses identified in the camel popula- 
tion within the B1 cluster representing contemporary drom- 
edary and human MERS-CoVs. The clustering with human 
MERS-CoVs known to have caused nosocomial outbreaks in 
the KSA, South Korea, and Jordan underlines the zoonotic 
potential of these camel-derived MERS-CoVs. 


Conclusions 


While the most recent nosocomial outbreaks in Jordan were 
linked to travel-related cases from KSA as reported by the 
WHO, the detection of B1-cluster-like MERS-CoV in drome- 
dary camels indicates that local introductions of MERS-CoV 
into the human population is a real possibility in Jordan, in 
addition to introductions via travel-associated cases. 
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